Spin-Hall edge spin polarization in a ballistic 2D electron system 
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Universal properties of spin-Hall effect in ballistic 2D electron systems are addressed. The net 
spin polarization across the edge of the conductor is second order, ~ A'^, in spin-orbit coupling 
constant independent of the form of the boundary potential, with the contributions of normal and 
evanescent modes each being ~ \/A but of opposite signs. This general result is confirmed by the 
analytical solution for a hard-wall boundary, which also yields the detailed distribution of the local 
spin polarization. The latter shows fast (Friedel) oscillations with the spin-orbit coupling entering 
via the period of slow beatings only. Long-wavelength contributions of evanescent and normal modes 
exactly cancel each other in the spectral distribution of the local spin density. 

PACS numbers: 73.23.-b, 72.25.-b 



Introduction. Spintronics addresses interplay of spin 
and orbital degrees of freedom in various transport, op- 
tical, etc. phenomena with the ultimate goal of achiev- 
ing spin manipulation in nanostructures. Special place 
in spintronics belongs to the spin-Hall effect predicted 
a long time ago [H, which recently entered the era of 
experimental observation 0, 0, Spin-Hall effect is 
characterized by a boundary (edge) spin polarization re- 
sulting when electric current is flowing through the sys- 
tem. It is customary classified into "extrinsic" (impurity- 
driven) SJE 0, Q and "intrinsic" (band-structure in- 
duced) [9|, llO| types. Initially theories of spin-Hall effect 
addressed such auxiliary quantity as spin current (for the 
review see Refs. ll[) in infinite systems, but later the 
emphasis shifted towards direct calculation of spin polar- 
ization in confined geometries. For diffusive systems the 
search is to complement the coupled spin-density diffu- 



sion equations 
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with suitable boundary conditions 



While it is now understood that in 2D systems spin- 
Hall effect generally occurs with more complicated spin- 
orbit couplings, any amount of disorder destroys spin- 
Hall effect in infinite systems with linear coupling 11|. 
It is, therefore, important to establish whether pure bal- 
listic systems (without disorder) can exhibit non-zero 
spin-Hall polarization. Driven by this motivation, studies 
of intrinsic spin-Hall effect in ballistic finite-size systems 
had been initiated, mostly by means of numerical meth- 
ods [23, im. It is significant to realize that the edge spin 

polarization in ballistic systems appears not as a result 
of electric field-driven acceleration of electron momenta 
(and associated with it precession of spins). As well 
known, electric field is absent inside an ideal ballistic con- 
ductor connected to reflectionless leads [22|. Spin-Hall 



spin accumulation in ballistic systems is due to the edge 
precession only. When the populations of left-moving 
and right moving states are different, the boundary scat- 
tering results in oscillatory (Friedel) edge polarization 
which is perpendicular to both the electric current and 
the normal direction to the boundary. Such polarization 



was considered numerically in Refs. [23| for a 2D elec- 
tron gas (2DEG). The case of a 3D hole semiconductor 
has also been analyzed recently A possibility of 

distinguishing edge effects from spin transport has been 
addressed experimentally in Ref. [25[. Edge spin polar- 



ization in parabolic quantum wires has been considered 
in Ref. |26i]. 

In the present paper we resolve analytically the bound- 
ary problem for a ballistic 2D electron gas with linear 
spin-orbit coupling ^27] and calculate the non-equilibrium 
edge spin-polarization in a wide strip connected to ideal 
leads with chemical potentials shifted by the applied volt- 
age. We present a general argument that the out-of plane 
spin polarization integrated over the lateral direction has 
a universal value, independent of the particular shape of 
the confining boundary potential U{x). In the limit of 
weak spin-orbit coupling, X <ti vp, 



Sz(x)dx 



X^eV 
12^2^' 



(1) 



where eV is the difference of the chemical potentials in 
the two leads, and vp the bulk value of the Fermi 
velocity (which is the same for both spin-split subbands). 

We then illustrate how this result arises from micro- 
scopic calculations in a model of a sharp boundary by 
obtaining the electron Green's functions in a concise an- 
alytical form. The obtained spin density is approximated 
by the expression (Ti = 1), 



eV 



cos {2mvpx) sin^ {mXx). 



(2) 



It is remarkable that the spin-orbit coupling constant en- 
ters via the period of beating only. 

Net spin polarization. Consider a semi-infinite ballistic 
2DEG described by the Hamiltonian 

&^ 

2m 



H 



dr -0^ 



iX{(7xdy 



d-ydx) 



U{x) 



V-, (3) 



where potential U{x) ensures boundary confinement (see 
Fig. 1). For the sake of simplicity we present derivation 
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U(x) 




2DEG 



FIG. 1: (Color online) Geometry of the system. Ideal leads 
filled by equilibrium electrons up to the chemical potentials 
shifted by the applied bias. The edge is formed by a confining 
potential U{x) vanishing for x — > oo. 



for the case of 'Rashba' spin-orbit interaction, though 
calculations for the 'Dresselhaus' coupling [2^ are com- 
pletely analogous |2^. The system is attached to two 
ideal reflectionless leads injecting equilibrium electrons 
into 2DEG. The chemical potentials of the leads are 
shifted by the applied voltage, eV. 

Since ky is an integral of motion (in case of reflec- 
tionless leads), it is convenient to use the Fourier rep- 
resentation along the y-axis for the electron operators, 
= Sfc Cky{x)e'^^^'^ . One can now derive the equa- 
tion of motion for the expectation value of the electron 
spin operator, s{ky,x) = i(cj. {x)acky{x)) , which can be 
readily written in the form, 



dtSy{ky,x) = -dj:Jl{ky,x) - 2\kySz{ky,x). (4) 

Here J| stands for the conventional operator of spin- 
current, i.e.. 



In a steady state the Ihs of Eq. (|4]) vanishes. Integrating 
Eq. (HI over the x-direction, we obtain for the net spin 
polarization. 



(5) 



It is straightforward to calculate the value of the {ky- 
resolved) spin current Jl^{ky, 00) inside the bulk of a 2D 
system: 



R Ul t. V / 



(6) 



/3=±1 k 



where np{kx,ky) stands for the population of different 
momentum states in the subband (3. Only "uncom- 
pensated" states contribute to the non-equilibrium spin 
polarization given by Eqs. ([5][6]); these states describe 
electrons that originate in the left lead {ky > 0) and 
belong to the energy interval near the Fermi energy, 
Ef < k'^/2m + I3k\ < Ep + eV . The integral ^ di- 
verges logarithmically at /cy — > 0. Assuming the same 
infrared cut-off in both subbands, k, we observe that the 




FIG. 2: (Color online) Dependence of the local spin polar- 
ization (|16() . in units of el/m/87r^, on the distance to the 
boundary for different values of spin-orbit coupling constant. 
Solid (red) line: \/vf = 0.1, dotted (blue) line: X/vf ~ 0.2, 
solid (black) line utilizes the approximate formula ([2} for 
X/vf ~ 0.2. The plot of Eq. (I17p is indistinguishable from 
the exact Eq. (|16|l on this scale. 



diverging In fc-contributions in the two subbands cancel 
each other, yielding in the linear (in V) response. 



Szdx 
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2A(27r)2 
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In 



(7) 



where vf — \J2Ef I'm, -|- is the Fermi velocity. Ex- 
panding this general result to the lowest non-vanishing 
order in XjvF we recover the net boundary polarization, 
Eq. (P). 

Evanescent modes . The reflection at the boundary 
mixes the two bulk subbands. Those states that belong 
to the domain, k~^ < ky < k~ , where k^ = m{vF T A), 
refereed to as evanescent states [i^ , are characterized by 
exponentially decaying contribution from the upper (-I-) 
subband. Repeating the calculations leading to Eq. ([7]) 
but now for the evanescent domain only, we obtain. 



"dx = 



eV 



2A(27r)2 




A _ In 1 + V^/'"F 

VF 1 - \J'X/VF 



(8) 



Remarkably, the net evanescent contribution is ^ a/A and 
is largely canceled by the contribution from the normal 
domain ky < k~^ , yielding Eq. ([1]) which is quadratic in A. 
This cancelation occurs for local spin density as well, see 
Eq. jll). 

Electron Green's function. Microscopic calculation of 
the local spin polarization can be most simply performed 
with the help of the electron Green's functions, 

six) = ^'^^T J ^[Gl^^ix,x)-GlFA^,x)]a, (9) 


where G§ ^{x,x') is the retarded Green's func- 
tion. Its advanced counterpart satisfies the condition 
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G^^^{x,x') = G^^^{x' ,x). The summation over energy 
in Eq. ([9]) is performed over the narrow strip of width 
eV, similar to Eqs. dUll). 

To illustrate how spin polarization arises from the so- 
lution of the Schrodinger equation, let us solve the prob- 
lem of a hard-wall boundary: U{x) = 0, for a; > and 
U{x) — oo ioT X < 0. The case of a smooth boundary 
where electrons adiabatically follow semiclassical trajec- 
tories fo r sp in-split subbands (30j will be considered sep- 
arately |31| . 

For a plane wave, ^ e**^"^, the equation for the Green's 
functions for {x,x' > 0) is 



2m 



\{i(Tydx + cF^ky) -I- E' 



G{x,x') 



-S{x 



(10) 

where the subscripts ky and E are omitted for simplicity, 
and E' = E ~ ky/2m. The boundary condition for the 
impenetrable wall is G{x,0) = G{0,x') = 0. We solve 
the problem by first noting that the following function 
C{x) satisfies both the homogeneous equation (flO)) and 
the boimdary condition £(0) = 0, 



Cix) 



1 



y-( 



(11) 



here * stands for the simple complex (not Hermitian) 
conjugate; the sum is taken over both subbands, with 
the projection matrix for the subband (3 defined as 



B 



/3 



1 



(12) 



where the absolute value of the electron momentum k^ 
is defined above Eq. ([5]) and its x-component is k^ = 

\ Kk^)"^ — fc^. Here we concentrate on the normal modes. 



where both k^ are real; rather simple modifications for 
the evanescent domain (where k^ is imaginary) are out- 
lined below. 

Using the function (fTT|) we can readily construct the 
solution for the inhomogeneous equation (fTO|) which sat- 
isfies the boundary condition G(0,x') — 0, 

G{x, x') = -2m[C{x)A{x') + 6(2; - x')C{x - x% (13) 

where A{x') is yet unknown matrix. Since Green's func- 
tion has to obey both the boundary condition G{x, 0) = 
and the equation conjugated to Eq. (jlOp . the matrix 
A{x') must be a homogeneous solution satisfying the con- 
dition A{Q) — —1. This determines it up to some con- 
stant matrix G different for the retarded and advanced 
Green's functions, i(x') = CC^x') - d^,t{x'), 

Gr,a{x,x') = -27n\c{x)CR^AC\x')-C{x)d^,t''(x') 

+Q{x - x')C{x- x')]. (14) 
The constant Gji {Ca) is most simply determined from 
the condition that the retarded (advanced) Green's func- 
tion does not contain the waves ~ e~^''^^ (e**^^^) propa- 
gating to (from) the boundary in the region x > x' . The 
calculations are straightforward but rather tedious. As a 
result one obtains, 



^a; ) T 2^ (k 



kx kx 



(15) 



where the upper (lower) sign corresponds to Gr{Ga)- 

Spin polarization. Making use of the derived Green's 
function we can now calculate the local spin polarization 
(HI). With the help of Eqs. (0 and ^ we obtain 



Szix) 



eV 



2{2T:)'^mVp [ Jq ky 



dk. 



^ (fc+Zc +kl- k+k^ ) si\i{2k+x) + sin(2fc^ x) - 2 sin((fc+ + k^)x) 



dky 

k+ ky 



fc^K(e*''- - e""^)^ - 2(fc+fc +ky + ik~ k) sin{k^ x){cos{k^ x) ~ e""^) 



(16) 



The first line here is the contribution of the normal modes 
while the second line comes from the evanescent modes, 

where k = \jky ~ {k+Y [13 • By calculating the integral 
over X it is straightforward to verify that the net contri- 
bution of the evanescent modes satisfies Eq. ([8]), being 
~ ^fX. This is mostly canceled by the contribution from 
the normal modes. The total net contribution of both 



the normal and evanescent states yields Eq. ([7]) in agree- 
ment with our general argument based on the equation 
of motion for spin operators, Eq. Behavior of the 
local spin density l|16p is shown in Fig. 2. In the most 
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relevant limit, X -4^ vp, Eq. ([T6| can be simplified to: 



Sz{x) = -eV 



mX 



VF Jo 



dy sin(a;yr^/Oe"''^^/« 



eV 



[sin {2k'^x) - 2k^x cos (2fc'^x)] 



ff^r /'^y sin[a;fc+(Vy+v/5>^)], (17) 



An^mVp jQ 

where ^"-"^ = 2m^/Xvp, and (5 = (fc^/fc+)^ — 1. Inte- 
grating this expression over a; we recover the net spin 
polarization ([T]). 

Spectral distribution of spin density. It is instructive 
to present the results in terms of the Fourier transform 
of the spin density, Sz{q) = 2 dxsz{x) sin qx. From 
Eq. p6|) we find Sz{q) in a form of piecewise continuous 
algebraic function defined in four domains. The surpris- 
ing feature of the spectral distribution revealed by this 
calculation is its vanishing, Sz{q) — 0, in the whole long- 
wavelength domain, < q < 2fc+. In particular, this 
shows the exact cancelation between normal and evanes- 
cent modes. For larger values of q we obtain to the lead- 
ing order in A, 



M = 



eVq 



-1, 2k+ < q < 2mvF, 
1, 2mvF <q<2k-, (18) 
-2/{qO\ 2k' < q. 



The plot of the spectral distribution is illustrated in 
Fig. 3. Remarkably, the net spin polarization (given by 
TT~^ J dqsz{q) / q) comes from the large-g tail (oc q'^) in 
the spectral density Sz{q). 

Conclusion . In this paper we solved analytically a 
problem of mesoscopic spin-Hall effect in a confined 2D 
electron system. We presented general arguments why 
the net spin polarization in a ballistic spin-Hall effect is 
independent of the boundary potential and confirmed the 
result by a straightforward calculation for the hard-wall 
boundary, for which the analytical solution was obtained. 
The spectral distribution of spin density consists of two 
narrow peaks of opposite sign whose heights are virtually 
independent of the small spin-orbit coupling constant. 
Surprisingly, long-wavelength contributions from evanes- 
cent and normal modes exactly cancel each other. Under- 
standing the level of universality of this result for arbi- 
trary boundary potentials remains a challenging problem. 
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supported by DOE, Award No. DE-FG02-06ER46313 
and by SFB TR 12. 
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